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1 Introduction 

Differential Galois theory and Lie symmetries are two different theoretical frameworks designed 
to deal with similar mathematical problems: the integration, reduction, classification and listing 
of solutions of differential equations. Both theories appeared simultaneously at the end of 19th 
century. However, the links between them remained hidden for a long time, mostly because of 
the apparent walls that separate mathematical disciplines. Differential Galois theory appears 
to be central to differential algebra. On the other hand, the theory of Lie symmetries belongs 
to the realm of local differential geometry. For a general exposition of both theories we refer 
the readers to [8, 21] and [11, 16] respectively. Throughout this paper, “differential Galois 
theory” refers to the Galois theory of systems of linear differential equations, also called the 
Picard-Vessiot theory. 

The works of Ziglin, Morales, Ramis, Sirno, Churchill and others show that differential Galois 
groups may measure obstructions to the existence of first integrals of Hamiltonian systems; 
regarding this approach to non-integrability of dynamical systems, see the recent survey [14] and 
references therein. This suggests an interplay between Lie symmetries and Galois group. Indeed, 
for Hamiltonian systems, the (musical) duality induced by the symplectic structure transforms 
a first integral into a Lie symmetry; so, obstructions to the existence of first integrals induce 
obstructions to the existence of Hamiltonian vector fields of Lie symmetries. 

In [2], Ayoul and Zung study a more general definition of integrability (Bogoyavlensky inte¬ 
grability [6] ), which includes directly Lie symmetries and they show again how the Galois groups 
of variational equations appear to give obstructions to the existence of such symmetries. 

*This paper is a contribution to the Special Issue on Algebraic Methods in Dynamical Systems. The full 
collection is available at http://www.emis.de/journals/SIGMA/AMDS2014.html 
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B. Malgrange has also suggested a link between Lie symmetries and his non-linear differential 
Galois theory (see, e.g., [12, Remark iii, p. 224]). The general idea is that differential Galois 
groups and Lie symmetries should be related in the sense that more symmetries should imply 
a smaller Galois group. Our results elaborate on this idea to turn it into precise statements 
which apply to systems of linear differential equations. 

In order to link Lie symmetries with differential Galois theory, we follow a geometrical ap¬ 
proach developed by the first two authors in [4], Some general results about symmetries were 
already stated in [5, Section 6], but in a more general context of automorphic systems. In par¬ 
ticular, it is implicit in [5] that the eigenring [3, 20] consists of vertical Lie symmetries. Our 
approach here is less abstract and more explicit (see Section 6.2). In connection with that, we 
point out that the relevance of the eigenring for the symmetries of linear differential equations 
was discovered independently by C. Jensen in the nice paper [9], but without any mention of a 
relationship with Picard-Vessiot theory. 

The results contained in this work may be summarized as follows: 

(a) The search for symmetries of systems of linear differential equations may be reduced to 
the search for a particular kind of symmetries, namely homogeneous polynomial vertical 
symmetries (Lemma 3.3, Proposition 3.5). 

(b) We focus on polynomial vertical infinitesimal symmetries and show that these are solutions 
of associated linear differential equations that can be deduced from the given equation (via 
the Tannakian theory). The differential Galois theories of those equations is related to 
that of our original system (Theorem 6.2). 

(c) The Galois group determines the Lie algebra of polynomial vertical symmetries with 
coefficients in the ground field (Theorem 6.4). 

(d) Each non-trivial vertical polynomial symmetry with coefficients in the ground field places 
constrains on the Galois group. Thus, the bigger the symmetry algebra, the smaller the 
Galois group. In several cases, a single symmetry may force the group to be abelian or 
solvable (Theorems 6.10 and 6.13, Corollaries 6.1 and 6.11). 

Section 2 contains the basic definitions of symmetries. Section 3 studies polynomial sym¬ 
metries and establishes part (a) above. Section 4 gives a dictionary between graded polynomial 
vector fields, linear actions and corresponding linear differential systems. A geometric definition 
of the differential Galois group is given in Section 5 and the comparisons with symmetries are 
derived in Section 6. 

The problem of linking differential Galois theory with Lie symmetries of differential operators 
has been studied by C. Athorne in [1] and by W.R. Oudshoorn and M. van der Put in [17]. This 
restriction to differential operators seems to complicate the relations between the Galois group 
and the Lie symmetries. So, in this work, we adopt a slightly different point of view. Instead 
of considering higher-order linear differential equations, we study systems of first-order linear 
differential equations. This leads us to a broader notion of infinitesimal symmetry which has 
an explicit relation with the differential Galois group. In Appendix A, we propose a comparison 
between the definitions of infinitesimal symmetries in the contexts of higher-order equations 
and first-order systems. 

2 Characteristic and vertical symmetries 

Let U be an open subset of the complex projective line CPi. By a function field K, we mean 
a subfield of the field of meromorphic functions on U such that /C contains the constants C and 
is closed with respect the derivation Clearly, fields of rational functions, elliptic functions, 
etc. are function fields. As shown by Seidenberg in [18, 19], any differential field which is finitely 
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generated over Q can be embedded into such a function field. We fix from now on the function 
field /C and its domain U. 

We consider a system of linear differential equations with coefficients in the function field K, 
^ = A i x )y, y = (yi,---,yn), A(x) e Qi(n,lC). (2.1) 

The poles of A(x) are called the singularities. The point at infinity may be considered as 
a singularity depending on its behavior after a suitable change of coordinates. The graphs of 
the solutions of ( 2 . 1 ) are the integral curves of the associated vector field: 

v d ^ . cl 

x = & + L <■« 

i,j =1 

which is a meromorphic vector field in U x C n . 

An infinitesimal symmetry of X is an analytic vector field Y defined in some open subset of 
U x C" such that [Y, X] and X are linearly dependent in their common domain of definition. 
In particular, vector fields of the form "/X where 7 is an analytic function defined in some open 
subset of IA x C n are infinitesimal symmetries of X. They are called characteristic symmetries 
of X. Since the definition of infinitesimal symmetries is local, we have sheaves of infinitesimal 
symmetries and characteristic symmetries of X in U x C n . 

The Lie bracket of two infinitesimal symmetries is also an infinitesimal symmetry. Hence, 
infinitesimal symmetries form a Lie algebra sheaf. Characteristic symmetries form an ideal of 
the Lie algebra sheaf of infinitesimal symmetries. An infinitesimal symmetry Y is a vertical 
symmetry if it is tangent to the fibers of the canonical projection IA x C n —> U , that is Yx = 0. 
Its expression in coordinates takes the form 
n fi 

r = Y < 2 ' 2 > 

If Y is a vertical symmetry, then the Lie bracket [X, Y] vanishes. 

The Lie algebra sheaf of vertical symmetries is canonically isomorphic to the quotient Lie 
algebra sheaf of all infinitesimal symmetries modulo the ideal of characteristic symmetries. If Y 
is an infinitesimal symmetry, we can take its vertical representative: Y = Y — (Yx)X. 

By this reduction, the algebra of vertical symmetries is isomorphic to the algebra of infini¬ 
tesimal symmetries modulo the ideal of characteristic symmetries. Thus, in order to study the 
symmetries of (2.1), it suffices to consider vertical symmetries. We consider symmetries that 
are defined in open subsets of the form V x C n with VC U. Such symmetries can be seen as 
sections of a sheaf defined in IA. 

Definition 2.1. The sheaf sym x in U assigns to each open subset VC U the Lie algebra of 
vertical infinitesimal symmetries of X defined in V x C n 

sym x {V) = {Y e £ an (V x C n ) \ \Y,X] = 0 and Yx = 0}, 
where £ an (V x C n ) stands for the Lie algebra of analytic vector fields in V x C n . 

Our objective is to describe the sections of this sheaf sym A - and its relation with the closed 
form solutions and the Picard-Vessiot theory of the system (2.1). From now on, when we mention 
a symmetry of X we mean a section of sym x , that is, a vertical infinitesimal symmetry. 

Our definition of infinitesimal symmetries is considered, for instance, in the reference 
book [11], where the vertical symmetries are called shuffling symmetries. We prefer our ter¬ 
minology, because it has a clearer geometrical meaning in our context of linear differential 
equations (i.e., fiber bundles, although we will not use this geometrical terminology explicitly 
here). Actually, the relation [X, Y\ = 0 sometimes appears in the literature as being the direct 
definition of an infinitesimal symmetry of X, see for instance [ 6 , 7]. 
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3 Polynomial vertical symmetries 


3.1 Polynomial vertical vector fields 

Definition 3.1. Let Y be a vertical vector field defined in V x C n with V an open subset of U 


Y = ^2fi(x,y) 

i —1 


d_ 

dyfi 


We say that Y is a polynomial vertical vector field when the components fi(x, y) are polynomials 
in the variables yi,... ,y n . 

Example 3.2. The (Euler) homogeneous vector field 




d 

dyi 


is polynomial vertical, indeed linear vertical. It is a symmetry of any system of linear differential 
equations. Hence, it is a global section of sym Y . 

The definitions of degree and homogeneous components of a polynomial vertical vector field 
are clear. Given a function field 1C of meromorphic functions on V C U, we can also speak of the 
polynomial vertical vector fields with coefficients in 1C. They are the polynomial vertical vector 
fields 


y = ^P i (x,y) 

i= 1 


d_ 
dyi ’ 


where Pi(x, y),..., P n (x, y) are in }C[yi, ...,y n ]- 


3.2 Homogeneous components of symmetries 

Let us consider a vertical infinitesimal symmetry Y £ syrn^V) with VC U. We can develop 
the Maclaurin series for the components of Y with respect the variables yi, ■ - ■ ,y n , obtaining 

analytic functions on V. We decompose Y as a sum of its 


polynomial vertical vector field of degree j (i.e., with respect 

with V CU be a symmetry of X. All the homogeneous compo¬ 
nents Yj of its Maclaurin series are symmetries of X. 

Proof. The Lie bracket can be computed componentwise because the map [X, •] is homogeneous 
of degree 0 in its action on vector fields, so we have 

0= \Y,X\ = [Y 0 ,X} + [Y ll X] + [Y 2 ,X} + -- - . 

For each j > 0, [A, Yj] is a homogeneous polynomial vertical vector field of degree j. Thus, all 
the terms of the above series vanish and we have proved the result. ■ 


Y = YI 9hoc{x)y c 


i= 1 ce£Z r } 


d_ 

dyfi 


where the functions gi, a {x) are 
homogeneous components 

y = y 0 + n + y 2 + • • •, 

where each Yj is a homogeneous 
to the y variables) in V x C”. 

Lemma 3.3. Let Y £ sym Y (V) 
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Remark 3.4. The Maclaurin series of Y depends only on the value of Y in small neighbourhood 
of V x {0} in V X C n . Assume that Y is a rational vertical symmetry of X, 


y = ^2fi{x,y) 

3 =1 


d_ 

dyi' 


fi{x,y) G C(x, yi ,.. .,y n ), 


whose polar set does not contain the curve CPi x { 0 } in CPj x C n (i.e., the denominators do 
not vanish indentically for y = 0). Then it admits a Maclaurin expansion in y and Lemma 3.3 
shows that each homogeneous component Yq, Yi, Y%, etc. of Y is a (homogeneous) polynomial 
vertical symmetry with coefficients in C(x). 


3.3 Homogeneous polynomial vertical symmetries 

The sheaf sym x contains the subsheaf of polynomial vertical symmetries that we denote by 
syrn^ 00 - Lemma 3.3 implies that the homogeneous components of polynomial vertical symmet¬ 
ries are also symmetries. Hence, we have a decomposition 

OO 

sym^°° = 0 symj, 

n=0 

where sym^ stands for the sheaf of homogeneous polynomial vertical symmetries of X of deg¬ 
ree r. These objects can be interpreted simultaneously in two complementary ways, as sheaves 
and as differential varieties: 

(a) As a sheaf, sym^ maps each open subset V C U to the set sym r Y (V) of homogeneous 
polynomial vertical symmetries of X defined in V x C n . 

(b) As a differential variety, sym Y maps each differential field extension /C C T to the set 
sym^(J r ) of homogeneous polynomial vertical symmetries of X with coefficients in T. 
Since the Lie bracket is computed algebraically, this makes perfect sense even if T is not 
a function field. Our forthcoming Theorem 6.1, stated and proved in Section 6 below, tells 
that symY is in fact a linear differential variety defined over 1C. 

If the function field /C contains the rational functions, then rational symmetries, as vector 
fields in CPi x C n , can be always reduced to polynomial vertical symmetries with coefficients 
in 1C: 

Proposition 3.5. Assume that 1C contains the field of rational functions C(x). Let Y be a ra¬ 
tional vector field in CPi x C n which is a non-characteristic rational symmetry of X and whose 
polar subset does not include the curve U x {0} inU x <C n . We consider the Maclaurin series 

Y - (Yx)X = T 0 + Yj + Y 2 + • • • , 

where each Y r is a homogeneous polynomial vertical vector field of degree r. Then, for each r, 
Y r G sym^/C) and for at least one index r, Y r is not zero. 

Proof. Lemma 3.3 and its Remark 3.4 show that the vector fields Y r are symmetries. We only 
need to check that they have coefficients in K. Let us consider the expression of Y in coordinates 



6 


B. Blazquez-Sanz, J.J. Morales-Ruiz and J.-A. Weil 


Then 


Y -{Yx)X = Y J H j {x,y) 

3 =1 



n 

Hj(x,y ) = ^2(fj(x,y) - h(x,y)a,ji(x)yi), 

i =1 


and 


n 


j =1 \a\=r 


d^Hj n .y a 8 
dy a ( X ’ K\d yj ' 


A direct examination of the expressions shows that they have coefficients in JC. ■ 

Since the Lie bracket is a graded operation, the sheaves sym^ are not in general Lie algebra 
sheaves. We have 


[ , ]: symj X sym^ ->• syrn^ 5 . 

For n > 1 only syrn^, sym^ and syrn^ © syniy are Lie algebra sheaves. Our next objective 
is to show that the sections of syrtij for each r are solutions of a hierarchy of linear differential 
systems canonically attached to (2.1). 


4 Polynomial vector fields in C n 

4.1 The Lie algebra of polynomial vector fields 

Polynomial vertical symmetries are polynomial vector fields along the fibers of the projection 
U X C n —> U. In this section we will give some remarks about the structure of the Lie alge¬ 
bra Jf[C n ] <0 ° of polynomial vector fields in C n 

X[C n ]<°° = j Pfiy ) A : Pfiy) € c[y 1} ... ,y n ] ) . 


» 2=1 


By taking homogeneous components, we have X[C n ] = 0^ o X r [C n ], where 

( n d 1 

X r [C n ] = < Pfiy) —— : Pfiy) homogeneous of degree r > . 

oyi 


. 2=1 


The Lie bracket respects the degree in the following way 
[,]: r[C n ]xl s [C"]4X r+s - 1 [C n ]. 

Remark 4.1. For n > 2, exactly two of the homogeneous components of X[C n ] are Lie subal¬ 
gebras: 

(a) The homogeneous component of degree zero X°[C n ]. It is the Lie algebra of the infinitesi¬ 
mal generators of the action of the group of translations in C n . It is an abelian Lie algebra 
canonically isomorphic to C n 


d 


£ U [C"]~C\ 

dy-j 


where {ei,..., e n } stands for the canonical basis of C n . 
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e eA y. 


£=0 


(b) The homogeneous component of degree one, X 1 [C n ]. It is the Lie algebra of linear vector 
fields in C n . It consists of the infinitesimal generators of the action of the group of linear 
transformations of C n . It is canonically isomorphic to g[(n, C) in the following sense: 

n f) 

X x [C n ] ~ fll(n,C), Y a b%'7T“ ^ A ’ 

i,j=i Vi 

where A stands for n x n matrix of entries a t j. 

Given an endomorphism A 6 gl(n, C), we let va denote its corresponding linear vector 
field in C n 

V"' <9 . . . \ d 

va := 2 ^ Wlth the identity va{v) = 

i,j=l 1 

We may easily check that this morphism is in fact an anti-isomorphism of Lie algebras: 
for any pair (A, B) of matrices, we have 

[v A ,V B \ +V[ A ,B] = 0. 

4.2 Induced linear actions 

Let us consider the canonical action of GL(n, C) on C n by linear transformations 
GL(n, C)xC n 4 C n , (A, y) ^ Ay , 

By abuse of notation, we denote non degenerate matrices and their associated linear transfor¬ 
mation on C n by the same symbols. Let A be a linear transformation and Y be a homogeneous 
polynomial vector field of degree r. 

Viewing A as a (linear) diffeomorphism of C n , we let A^(Y) denote the vector field trans¬ 
formed by A. In general, for any diffeomorphism F, we let F*(Y)(F(p)) = dF(Y(p)), i.e., 

F*(Y)(p) = dF(Y(F- 1 (p))). 

This defines a natural action of diffeomorphisms of C n on vector fields of C n . 

It is easy to check that A*(Y) is also a homogeneous polynomial vector field of degree r. 
Thus, for each r > 0 we have an induced representation 

T r : GL(n,C) —► GL(3T[C n ]), $ r (A)(Y) = A*(Y), 

which yields a linear representations of GL(n, C) on the finite-dimensional vector spaces X r [C n ]. 
This action can be differentiated at the identity obtaining the infinitesimal action 

d 


K- 0[(n,C)^End(3r[C n ]), $J.(A)(Y) = 


de 


6=0 


T(e £ ^)(V). 


The following remark is key to connect the definition of symmetry with the differential Galois 
theoretic aspects of equation (2.1). 

Lemma 4.2. The infinitesimal action T 7 o/gl(n, C) in 3f r '[C Tl ] coincides up to a change of sign, 
by the canonical isomorphism between X x [C n ] and gl(n,C), with the Lie bracket action of linear 
vector fields on X r [C n ] . That is, for any endomorphism A and homogeneous polynomial vector 
field Y in C n 

<P' r (A)(Y) =-[v a ,Y]. 

Proof. Let us define, for each A and e, o £ : C n —> C n , the map that sends each y € C n to e eA y. 
Thus, {cr £ } £e c is the flow of the vector field va- We have 


<*Mm = f £ 


£=0 


^ a )(Y) - I 


<?e*(Y) = -Lie^y = —[va, L"], 


e=0 


by the usual geometric definition of Lie derivative. 
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5 Lie—Vessiot hierarchy and Galois group 

5.1 The Lie—Vessiot hierarchy 

Here we recall some of the definitions from [5], adapted to the particular case of linear equa¬ 
tions. Let E be a finite-dimensional complex vector space and (E, 'I') a linear representation 
of GL(n, C) in E. The group morphism T induces a Lie algebra morphism 'E: 

T: GL(n,C) GL(£), V: 0 l(n,C) -> End(£). 

This morphism transports our linear differential system (2.1) to a linear differential system in 
U x E with coefficients in 1C: 


( 77? 

- = T'(A(x))(u). 


(5.1) 


We say that system (5.1) is the Lie-Vessiot system induced by (2.1) in the representation {E, 'L). 
They are the geometric analog of the differential systems obtained by Tannakian correspondence 
on tensor constructions in standard differential Galois theory, see [10, 21]. 

A solution of the Lie-Vessiot system (5.1) in V C U is an analytic map V —> E that satisfies 
the equations. Given a differential field extension 1C C E, a solution of the Lie-Vessiot sys¬ 
tem (5.1) in T can be thought of as being an element of E < 8 >c J~- If we take a basis v\,... ,v n 
of E and denote by Ai,...,A n their corresponding linear coordinate functions, the functions 
bij(x ) = Aj(T , (A(x))uj) are elements of 1C and the differential equation can be written in coor¬ 
dinates 


dX 

dx 


= B(x) A, A = (Ai,..., A m ), B(x) G 0 [(m,/C). 


(5.2) 


are 


Note that, as the Lie-Vessiot construction is a Lie algebra morphism, the poles of B{x 
exactly the poles of A(x). 

There is a natural relation between solutions of (2.1) and of its induced Lie-Vessiot sys¬ 
tem (5.1): if M(x) is a fundamental matrix of solutions of (2.1) defined in V C U then, for all 
vq G E, v(x) = '&(M(x))v o is a particular solution of the Lie-Vessiot system (5.1). 


5.2 The Galois group 

The Lie-Vessiot systems induced by (2.1) form a hierarchy of differential equations which encode 
the differential algebraic properties of (2.1). It allows us to give a “geometric” definition of the 
differential Galois group. In this definition we are concerned with two kind of solutions of the 
Lie-Vessiot systems: 

(a) We say that a solution v(x) of the Lie-Vessiot system (5.1) is IC-rational if it belongs to 
E <8>c V. This means that v(x) has its coordinates (5.2) in 1C. 

(b) We say that an element w(x) of E®clC is a 1C-exponential pre-solution of the Lie-Vessiot 
system (5.1) if there is a function b(x) G 1C such that 

^ — T / (A(x))(rc) = —b(x)w. 
dx 

It models the case in which v(x) = exp (f b(x)dx)w(x ) is a solution of the Lie-Vessiot 
system (5.1) or, equivalently, the class (v(x)} is a rational solution of the equation in the 
projective space P(i?) obtained by reducing the Lie-Vessiot system (5.1) by the Euler 
homogeneous vector field (see Example 3.2) of symmetries. 
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Note that the concept of /C-exponential pre-solution extends that of /C-rational solution: any 
/C-rational solution is a /C-exponential pre-solution in which the multiplier b(x) vanishes. 

Definition 5.1. Let us fix an xo £ U which is not a singularity of (2.1). We say that a non¬ 
degenerate matrix o E GL(n, C) is Galoisian at xq if for any linear representation it 

satisfies the two following conditions: 

(a) For any /C-rational solution v(x) of any induced Lie-Vessiot system (5.1), 'F(cr)(u(xo)) = 
v(xq). Note that, if v(x) is a /C-rational solution then x(xo) is well defined as an element 
of E. 

(b) For any /C-exponential pre-solution w(x) of any induced Lie-Vessiot system (5.1), for 
which w(xo) is well defined, w(xo) is an eigenvector of 'F(er). 

The Galoisian matrices at xq form a group Gal(xo, A), called the Galois group of (2.1) at the 
point xo- It is the stabilizer of all the values at xo of /C-rational solutions, and the lines spanned 
by the values at xq of /C-exponential pre-solutions of induced Lie-Vessiot systems. 

Although this geometric definition may seem different from the standard ones from Picard- 
Vessiot theory, it produces the same group. Choose a normalized local solution matrix at xo, i.e., 
one with initial condition being the identity at the point xo; then v{xq) will be the coordinates 
of the invariant v(x) on this normalized basis of solutions. Our definition hence says that a mat¬ 
rix a is in Gal(xo, A) if and only if it admits all the (semi-)invariants of the (Picard-Vessiot) 
differential Galois group as (semi-)invariants. 

The following facts are well known in differential Galois theory (we refer the interested 
reader to [8, 21] for a general exposition, or to [4, 5] for an exposition which is consistent with 
our geometric definition): 

(a) The Galois group Gal(xo, A) is an algebraic subgroup of GL(n, C). 

(b) The Galois groups at two different non-singular points xq and x\ are conjugated. We will 
write Gal (A) to denote this abstract Galois group, that we call the Galois group of the 
equation (2.1) over /C. 

(c) The system (2.1) is integrable by Liouvillian functions if and only if Gal(xo, A) is a virtually 
solvable group, i.e., its identity component is solvable. 

The next lemma encodes the expected Galois correspondence. The reader may check [4, 
Proposition 5.4] for a geometric proof that relies on Lie’s reduction method and Chevalley 
theorem. 

Lemma 5.2. Let (£\\k) be a linear representation o/GL(n, C) and let z{x) be a solution of its 
corresponding induced Lie-Vessiot system (5.1). Then, z(x) is a JC-rational solution if and only 
if for all Galoisian matrices at xo, <x £ Gal(xo, A), we have Tf(cr)z(x o) = z(x o). 

6 Symmetries vs Galois 

6.1 Symmetries as solutions of the Lie—Vessiot hierarchy 

A key of the relation between the Galois group and the symmetries is the fact that polynomial 
vertical vector fields in V C U of can be seen as maps V —> X[C n ]. For a given polynomial vertical 
vector field Y and xo 6 V we will write Y (xo) for the value ofY at xo- It is a polynomial vector 
field in C n that corresponds to the restriction of Y to the fibre {xo} x C n . It is clear that, 
for general vertical vector fields, [Y,Z](xo) = [V(xo), Z{xq)\. If we restrict our considerations 
to homogeneous polynomial vertical vector fields of a fixed degree r, then X r [C n ] turns out 
to be a finite-dimensional complex space. This will allow us to describe polynomial vertical 
symmetries as solutions of some systems of the Lie-Vessiot hierarchy. 
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Theorem 6.1. Let Y be a homogeneous polynomial vertical vector field of degree r defined 
in V x C n with VC U. Then, Y is a symmetry of (2.1) if and only if, as a map from V 
to X r [ C n ], it is a solution of the Lie-Vessiot system induced in the representation (<F r , £ r [C"]) 
from Section 4.2. 

In other words, Y £ sym^(V) if and only if it satisfies 


dY 

dx 


K(A(x))Y. 


( 6 . 1 ) 


Proof. Recall that if va = va{x) is the linear vertical vector field in U X C n corresponding to 
the matrix A(x), then X = + va(x). Let us compute the Lie bracket 


[X, Y} 


d_ 

dx 


+ v a (x),Y 


dY 

dx 


+ [va(x),Y]. 


Thus, Y is a symmetry if and only if '-jA = — [va(x), Y]. Finally, by Lemma 4.2 we have that Y 
satisfies the stated differential equation if and only if it is a symmetry. ■ 


If we denote by X[C n ] <0 ° the polynomial vector fields in C n , then: 

Corollary 6.2. Let xq be a non-singular point, and V a simply-connected neighbourhood of xq 
in U. Then, for each polynomial vector field W 0 ' £ X[C n ] <0 ° there is a unique polynomial ver¬ 
tical symmetry Y £ sym x (V) <0 ° such thatY(x o) = Moreover, sym x (V) <0 ° andX[ C n ] <oc 

are isomorphic Lie algebras. 

Proof. The map symjjf(V) <0 ° —> X[C n ] <0 °, Y t —> Y(x o), is a Lie algebra morphism since 
the computation of the Lie bracket and the restriction to the fiber {xo} x C n are commuting 
processes. We have to see that it is an isomorphism. Let r be the degree of and 

y(°) = y Q (0) + yW + . . . + yW 

the decomposition of Y in homogeneous components. Let Y^ be the solution in V of the 
Cauchy problem 

^ = & j (A(x))Y k , Y k (x 0 ) = Y^. 

The existence and uniqueness of the solution guarantees that 


Y = T 0 + Yi + ■ ■ ■ + Y r 

is the only polynomial vertical symmetry such that Y(xo) = V®. ■ 

We will now write the system (6.1) in more explicit form using tensor products. 

Proposition 6.3. Let N := ( ra+ ™'~ 1 ) denote the number of monomials of degree rn in n va¬ 
riables. We define the matrix A rn ■= A <S> Id a? + Id„, <g)sym''' , '(A*) of size nN. 

The system y' = A m y has a rational solution 

y = [«i,i(x),... ,ajv,i(x),... ,oij(x),.. .,a N}j (x ),..., 

,n{x), ■ • ■ , <2AT,ra(x)] , Ojj(x) £ /C 

if and only if X admits the homogeneous vertical symmetry 
n / N 

Y = ^2 ( ^2 a id (x)m(yi,..., y n ) 

j =i \i =l 

where Hi{y\, ■ ■ ■, y n ) denotes the i-th monomial (for the lexicographic order) of degree m in the n 
variables yi ,..., y n . 
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Proof. Let ( Ai,d ) denote the differential module associated with y' = Ay. Letting yi,... ,y n 
denote a basis of the dual M*, we see that X represents the action of d on A4* ( X acts on 
the first integrals rather than on the solutions, see [13, 15, 22]). The monomials of degree rri in 
the yi,... ,y n form a basis of Sym m (Ad*). So, a homogeneous vertical polynomial symmetry is 
a map from A4* to Sym m (Af*) which furthermore commutes with X (and hence with d ). Now 

Horn (. M *, Syrn m (M*)) = (M*)* <8> Syrn m (M*) = M <8> Syrn m (M*), 

so that 

Homg (M*, Sy = ker (d , M ® Sym m (M*)). 

This shows that the coefficients of Y are exactly the rational solutions of y' = Amy. ■ 

6.2 Symmetries vs Galois 

The intrinsic relation between the Galois group and the Lie algebra of symmetries of (5.2) is 
made explicit by the following result. 

Theorem 6.4. Let be a polynomial vector field in C n , and xo a non-singular point of (2.1). 
There is a polynomial vertical symmetry Y of X with coefficients in K, such that Y(x o) = 
if and only if for each Galoisian matrix a £ Gal(xo,^C), ct^YA) = y(°). 

Proof. It follows from our definition of Galois group. By Theorem 6.2, for each one of them 
there is a polynomial vertical Y symmetry such that Y(x o) = defined in a neigbourhood 
of xo. By Theorem 6.1 this is a solution of a Lie-Vessiot system induced by X, here we consider 
all the homogeneous components simultaneously. Finally, by Lemma 5.2, this solution has 
coefficients in /C if and only if A*Y(xo) = Y(x o) for all Galoisian matrices at xq. ■ 

Hence, the Galois group Gal(xo,^0 of (2.1) determines the Lie algebra sym <0O (/C) of poly¬ 
nomial vertical /C-rational symmetries in the following sense. The Lie algebra sym <0O (/C) is 
isomorphic to X[C n ] Gal ^°’ A ^, the Lie algebra of polynomial vector fields fixed by the action of 
Gal(xo,^0 by linear transformations in C n . However, we do not have a reciprocal: in general, 
the Galois group is contained in the stabilizer of the Lie algebra of /C-rational symmetries. 

It is possible also to dualize the situation and to consider the Galois group itself as symmetries 
of the infinitesimal symmetries of equation (2.1) as follows. Given a polynomial vector field Y 
in C n a linear symmetry of Y is a non-degenerated matrix a such that afiY = Y. Here, a stands 
for the transformation 

If the expression in coordinates of Y is 

then the matrix cr is linear symmetry of Y if and only if it satisfies the equations 

n 

Piipy) = ( 6 - 2 ) 

3 = 1 

If we are looking for the linear symmetries of a homogeneous polynomial vector field of degree m, 
it yields a total of n x ( n+ ™~ ) equations. Thus, for generic polynomial vector fields of high 
degree the group of linear symmetries reduces to the identity. 


y = (2/1, ■ ,y n ) -> <^y = 


vijVj’ ■ ■ ■ > 


3 = 1 
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Example 6.5. Let us compute the symmetries of the quadratic vector field 
Y-v 2 — 

in C 2 . Equations (6.2) for this particular case yield 


oil 2/i + 20210222/12/2 + crlivl = <7ll2/|, 


0 = 0 - 212 / 2 , 


equating each coefficient, we obtain 

011 = 0I2, 021 = 0. 


Thus, the group of linear symmetries is 




A e C\ 



(6.3) 


The Lie-Vessiot induced system for polynomial vertical symmetries of arbitrary degree r is 


dY 

dx 


~[va(x),Y]. 


We can consider all the homogeneous components simultaneously. Theorem 6.4 can be restated 
in the following terms: 

Corollary 6.6. Let Y^ be a polynomial vector field in C n , and xq a non-singular point of 
equation (2.1). The necessary and sufficient condition for the existence of a polynomial vertical 
IC-rational symmetry Y of X such that Y(x o) = is that the Galois group Gal(xo,^0 is 
contained in the group of linear symmetries of Yq • 


Example 6.7. Let us consider the system 


dyi 

dx 


2a(x)yi + b{x)y 2 , 


dy 2 
dx 


a(x)y 2 , 


where a(x), b(x), are arbitrary functions in /C. A direct computation of the Lie bracket says that 
Y = y{fiy 2 is a symmetry, and thus the Galois group of the equation (for any function field 1C) 
is contained in the group (6.3). 


Let us check now how the polynomial vertical symmetries of degrees zero and one look like, 
and what kind information about the Galois group they carry. 


6.2.1 Symmetries of degree zero 

The canonical isomorphism X°[C n ] ~ C n stated in Section 4.1 tell us that the linear representa¬ 
tion (4>o, Jf°[C Tl ]) is just an isomorphism. In particular, if 

n P) 

y = 

ST d «‘ 

then 

n rl 

&(A(x))Y = a i j (x)f j {x) — , 
i,j =1 Vl 

and thus we have the following result: 
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Proposition 6 . 8 . A n-tuple of functions y = (</> i(x),..., cj> n (x)) is a solution of ( 2 . 1 ) if and 

n 

only ifY is a symmetry of (2.1). 

i =1 

This proposition can be understood as an infinitesimal version of the superposition principle. 
If y(x) and (j>(x) are solutions of (2.1), then for all e, y(x) + £(j)(x) is also a solution. For 
a fixed cj)(x) and a general y(x), the derivative of this monoparametric family of solutions with 
respect to e is a vertical vector field, namely, the symmetry Y. 

Remark 6.9. Given a symmetry of degree zero of (2.1) (i.e., a solution), we can reduce the 
system of n differential equations ( 2 . 1 ) to a system of n — 1 differential equations, by means of 
a suitable gauge transformation. This simple observation can be viewed as a generalization of 
the classical result of d’Alembert: the order of a linear differential equation can be reduced by 
one when a particular solution is known. 

6.2.2 Linear symmetries 

Homogeneous polynomial vertical symmetries of degree one are called linear symmetries. The 

_» 11 

homogeneous vector field h = and its multiples gives us a trivial monoparametric family 

1=1 

linear vertical symmetries for any system of differential equations. The canonical isomorphism 
£ 1 [C Tl ] ~ g[(n,C) stated in Section 4.1 tell us the that Lie-Vessiot system induced by the 
representation (<J>i, £ 1 [C n ]) can be seen as a matrix equation. If we write 

n d 

where B = (bij ) stands for a n x n matrix of undetermined functions, the induced system is 
written as 

f =[«),«]■ (6.4) 

This is the equation of isospectral deformations induced by A{x) and has been exhaustively 
studied. The set of rational solutions of (6.4) is called the eigenring , see [3] for an extensive 
study of its properties, notably to decompose linear differential systems. If B(x) is a solution 
of (6.4) it is well known that the Jordan canonical form of B(x ) does not depend on the point x. 
Thus, given a linear symmetry Y with matrix B(x), we will classify it according to its Jordan 
canonical form: 

(a) If B(x) has at least two different eigenvalues we will say that Y is a decomposer symmetry. 

(b) If all the eigenvalues of B{x) are different we will say that Y is a complete decomposer 
symmetry. 

(b) If the eigenspaces of B(x) are one-dimensional we say that Y is a solver symmetry. That 
means that its Jordan canonical form does not contain any block of the form 

A 0\ 

0 XJ ' 

The following theorem is very close to some of the results of C. Jensen in [9, Section 9], on 
integration by quadratures, and results of M.A. Barkatou in [3] on decomposition - although we 
obtain it by different means and relate it with the linear symetries and the Galois group of the 
system. 
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Theorem 6.10. The following are equivalent: 

(а) There is a decomposer symmetry Y E sym(/C) with k different eigenvalues of multiplicity 

n,...,r k . 

(б) The Galois group Gal(xo,A!) is conjugated to a sugbroup of the group of block-diagonal 
matrices GL(ri,C) x ••• x GL(r*;,C). 

Proof, (a) ==> (b). Let .To be a non-singular point and B the matrix defined by Y(xq), with 
eigenvalues X\,... ,Xk of multiplicities r\,... ,rk- Let us consider the decomposition 

C n = Ei © • • • © E k , (6.5) 

where the spaces Ej = ker(E — Aj Id) are the generalized eigenspaces of B. The group 

G = {<tG GL(n, C): a{Ei) = Ei for all* = 1,..., k} 

is clearly conjugated to the group of block-diagonal matrices. Let us see that all Galoisian 
matrices at to are in G. If a is Galoisian then, ct*(T(to)) = Y(to), but that means aBa~ 1 = B , 
so a conjugates B with itself, and thus it sends generalized eigenspaces of B to themselves. 

(b) =>■ (a). Let us assume that Gal(To, X) is conjugated to a subgroup of the group of block- 
diagonal matrices. Then, we have a decomposition of C n in subspaces as in formula (6.5), such 
that for all a E Gal(To, X), and index i = 1,..., k, a(Ej.) = E). Let us consider the following 
linear vector fields hi in C n for i = 1 ,..., k defined by properties: 

hi | = h, where h stands for the Euler homogeneous vector field, 

h i \E j =0, tfi^j- 

Let us consider m ,..., fik different complex numbers and define 

Y {0) = ^2 mhi. 

1=1 

is stabilized by any Galoisian matrix and then, by Lemma 5.2, there is /C-rational sym¬ 
metry Y whose value at to is Y - ®. This symmetry Y is the decomposer symmetry of the 
statement. ■ 

Corollary 6.11. The following are equivalent: 

(a) There is a complete-decomposer linear symmetry in sym^(/C). 

( b ) The Galois group Gal(To, X) is conjugated to subgroup of the group of diagonal matrices 
(( C*) n C GL(n,C). 

Proof. The statement is the particular case of Theorem 6.10 in which all the eigenvalues are 
simple. ■ 

Remark 6.12. The existence of a complete-decomposer linear symmetry implies the existence 
of a n-dimensional abelian Lie algebra of symmetries. Let us consider a complete-decomposer 
linear symmetry Y and To a non-singular point. Let B be the matrix of Y(to), and vi,...,v n 
be a basis of eigenvector of B. As before, we define vector fields: 

hi{vf) = h[vi ), where h stands for the homogeneous vector field, 

hi(vj) = 0, if i^j. 

It is easy to check that the matrices of the vector fields hi have common eigenvectors and then 
[hi,hj] = 0. For all Galoisian matrix a at To we have cr*(hi) = hi, and thus by Lemma 5.2 
there are linear vertical /C-rational symmetries H ±,..., H n such that Hi{ to) = hi. They form 
a n-dimensional abelian Lie algebra. 
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Theorem 6.13. If there is a solver symmetry Y € sym(/C) then the Galois group Gal(:ro,Al) 
is conjugated to a subgroup of the triangular group. 

Proof. Let us consider first the case in which B, the matrix of Y(xq), has only one eigenvalue A 
of multiplicity n. For each Galoisian matrix a at xq we have aBa~ l = B. We have a chain of 
subspaces 

0 C ker (B - Aid) C ker (B - Aid) 2 C • • • C ker (B - Aid)”" 1 C C n . 

In general, a (ker (B — Aid)- 5 ) = ker(crL>er _1 — Aid)- 7 and thus Galoisian matrices respect the 
chain of subspaces. In other words, they are triangular matrices in some suitable basis. For 
the general case, with different eigenvalues, we first consider the decomposition of the group by 
block-diagonal matrices given in Theorem 6.10, and then we apply the above argument. ■ 

The results of this subsection supported Lie’s idea that Lie symmetries are useful for the 
integrability of the differential equation (2.1) by quadratures or at least its partial integrability 
or reduction: 

(1) By Remark 6.9, the existence of a symmetry of order zero implies a reduction of the order. 

(2) The decomposition in block-diagonal form of the Galois group implies that, by means of 
a gauge transformation, we can transform the equation (2.1) in a direct sum of linear 
differential equations (Kolchin or Lie-Kolchin reduction, see [5, 21]). By Theorem 6.10, 
we fall in this case for a decomposer symmetry. 

(3) By Remark 6.12 and Theorem 6.13, the existence of either a complete-decomposer or 
a solver symmetry implies the solvability of the equation by Liouvillian functions. 

Also it is not difficult to obtain some results for the Hamiltonian symmetries in the symplectic 
case, i.e., for non-autonomous linear Hamiltonian systems. In some sense, this approach would 
shed light on the references [13, 15] from the Lie point of view. 


A Symmetries of higher-order equations vs first-order systems 


There are two different ways to present the theory of linear differential equations. The first one 
deals with a single higher-order linear differential equation: 


d n y 


d^y 


H-b a 0 (x)y = 0. 


The second one deals with a system of first-order linear differential equations 



/ 2/o ^ 


/ 

0 

1 

0 ... 

0 \ 

d 

2/1 



0 

0 

1 ... 

0 

dx 









\yn-i) \-a 0 {x) -ai(x) . -a n - 1 (x)J 


^ 2/o ^ 
2/1 

\yn-l) 


(A.l) 


(A.2) 


The system (A.2) is called the companion system of (A.l). The variable yi represents the 
i-th derivative of the function y. In this paper, we have considered systems instead of higher- 
order equations. There is a pragmatic reason: although the differential Galois theories for 
equations (A.l) and (A.2) are the same (see, for instance [21, Section 2.1]), their Lie symmet¬ 
ry theories are not. Higher-order equations have less symmetries than systems of first-order 
differential equations. 










16 


B. Blazquez-Sanz, J.J. Morales-Ruiz and J.-A. Weil 


By definition, an infinitesimal point symmetry of (A.l) is a vector field 

_ _ . / , d d 

Y = ^v)g x +v(*,y)g y 

in the plane x, y whose flow maps solutions to solutions. Let us consider the vector field 

d d d 

X = ~dx +Vl lkjo^ -( a o(x)yo H- 1 - a n -i(x)y n -\) —— 

The vector field Y extends to a unique vector field Y in the jet space of coordinates x,yo, 

■ ■ ■ ■> Un-i 

0 0 Q 0 

y = ^(x,y)— +T]{x,y 0 )~ -hr?i(x,2/o,yi)w H-h Vn-i(x, y 0 ,..., y n -i) w-, 

ox oyo o oy n -1 


satisfying the conditions (see [16, Section 2.3]) 


[y,x] e (x), 

Uey(dyo - yidx,..., dy n - 2 - y n -idx) C ( dy 0 - y±dx ,..., dy n - 2 - y n -idx). 

It is known that the Lie algebra of point symmetries (in some open subset) of a linear differential 
equation of order > 2 is finite-dimensional. The above conditions allow us to generalize the idea 
of infinitesimal point symmetry. Any vector field in the jet space of coordinates x,yo,, y n -i 
is called an infinitesimal contact symmetry if is satisfies 


[Z,X\ € (X), 

Lie z (dy 0 - y\dx ,..., dy n - 2 - y n -idx) C ( dy 0 - yidx, ..., dy n - 2 - y n -idx). 

Infinitesimal point symmetries form a Lie subalgebra of the Lie algebra of infinitesimal contact 
symmetries. 

On the other hand, an infinitesimal symmetry of the system (A.2) is a vector field Z such 
that [Z,X\ G (X). See for instance [1] and [11, pp. 12-16]. It is clear that the Lie algebra 
of symmetries of the companion system (A.2) contains the Lie algebra of infinitesimal contact 
symmetries of the higher-order differential equation (A.l). 

In this paper, we have explored the relation between some Lie algebras of symmetries of a first- 
order system (A.2) and its differential Galois group. The relation between the Lie algebra of 
infinitesimal point symmetries of an operator (A.l) and its differential Galois group has been 
investigated, with rather negative results, in [1] and [17]. 
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